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Abstract 



We consider theories with gravity, gauge fields and scalars in four- dimensional 
asymptotically flat space-time. By studying the equations of motion directly 
we show that the attractor mechanism can work for non-supersymmetric 
extremal black holes. Two conditions are sufficient for this, they are conve- 
niently stated in terms of an effective potential involving the scalars and the 
charges carried by the black hole. Our analysis applies to black holes in the- 
ories with Af ^ 1 supersymmetry, as well as non-supersymmetric black holes 
in theories with M = 2 supersymmetry. Similar results are also obtained 
for extremal black holes in asymptotically Anti-de Sitter space and in higher 
dimensions. 
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1 Introduction 

Black holes in M = 2 supersymmetric theories are known to exhibit a fas- 
cinating phenomenon called the the attractor mechanism. There is a family 
of black hole solutions in these theories which are spherically symmetric, 
extremal black holes, with double-zero horizons l . In these solutions sev- 
eral moduli fields are drawn to fixed values at the horizon of the black hole 
regardless of the values they take at asymptotic infinity. The fixed values 
are determined entirely by the charges carried by the black hole. This phe- 
nomenon was first discussed by [T] and has been studied quite extensively 
since then OE1I110EIIZ1IHIE1 UD] • It has gained considerable attention re- 
cently due to the conjecture of [TTj and related developments [T2 l 1X3 } IT U IT5] . 

So far the attractor phenomenon has been studied almost exclusively in 
the context of BPS black holes in the M = 2 theories. The aim of this paper 
is to examine if it is more general and can happen for non-supersymmetric 
black holes as well. These black holes might be solutions in theories which 

1 By a double-zero horizon we mean a horizon for which the surface gravity vanishes 
because the goo component of the metric has a double-zero (in appropriate coordinates), 
as in an extremal Reissner Nordstrom black hole. 
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have no supersymmetry or might be non-supersymmetric solutions in M ^ 1 
super symmetric theories. 

There are two motivations for this investigation. First, a non-supersymmetric 
attractor mechanism might help in the study of non-supersymmetric black 
holes, especially their entropy. Second, given interesting parallels between 
flux compactifications and the attractor mechanisms, a non-supersymmetric 
attractor phenomenon might lead to useful lessons for non-supersymmetric 
flux compactifications. For example, it could help in finding dual descrip- 
tions of such compactifications. This might help to single out vacua with a 
small cosmo logical constant. Or it might suggest ways to weight vacua with 
small cosmological constants preferentially while summing over all of them 
2 . These lessons would be helpful in light of the vast number of vacua that 
have been recently uncovered in string theory [To] . 

An intuitive argument for the attractor mechanism is as follows. One 
expects that the total number of microstates corresponding to an extremal 
black hole is determined by the quantised charges it carries, and therefore 
does not vary continuously. If the counting of microstates agrees with the 
Bekenstein-Hawking entropy, that is the horizon area, it too should be de- 
termined by the charges alone. This suggests that the moduli fields which 
determine the horizon area take fixed values at the horizon, and these fixed 
values depend only on the charges, independent of the asymptotic values 
for the moduli. While this argument is only suggestive what is notable for 
the present discussion is that it does not rely on supersymmetry. This pro- 
vides further motivation to search for a non-supersymmetric version of the 
attractor mechanism. 

The theories we consider in this paper consist of gravity, gauge fields and 
scalar fields. The scalars determine the gauge couplings and there by couple 
to the gauge fields. It is important that the scalars do not have a potential 
of their own that gives them in particular a mass. Such a potential would 
mean that the scalars are no longer moduli. 

We first study black holes in asymptotically flat four dimensions. Our 
main result is to show that the attractor mechanism works quite generally 
in such theories provided two conditions are met. These conditions are suc- 
cinctly stated in terms of an "effective potential" V e ff for the scalar fields, 
(pi. The effective potential is proportional to the energy density in the elec- 
tromagnetic field and arises after solving for the gauge fields in terms of the 
charges carried by the black hole, as we explain in more detail below. The 
two conditions that need to be met are the following. First, as a function 

2 For a recent attempt along these lines where supersymmetric compactifications have 
been considered, see ■ 



3 



of the moduli fields V e ff must have a critical point, diV e ff((f>io) = 0. And 
second, the matrix of second derivatives of the effective potential at the criti- 
cal point, dijV e ff(4>io), must be have only positive eigenvalues. The resulting 
attractor values for the moduli are the critical values, (p i0 . And the entropy 
of the black hole is proportional to V((f> i0 ), and is thus independent of the 
asymptotic values for the moduli. It is worth noting that the two conditions 
stated above are met by BPS black hole attractors in an M = 2 theory. 

The analysis for BPS attractors simplifies greatly due to the use of the 
first order equations of motion. In the non-supersymmetric context one has 
to work with the second order equations directly and this complicates the 
analysis. We find evidence for the attractor mechanism in three different 
ways. First, we analyse the equations using perturbation theory. The starting 
point is a black hole solution, where the asymptotic values for the moduli 
equal their critical values. This gives rise to an extremal Reissner Nordstrom 
black hole. By varying the asymptotic values a little at infinity one can 
now study the resulting equations in perturbation theory. Even though the 
equations are second order, in perturbation theory they are linear, and this 
makes them tractable. The analysis can be carried out quite generally for any 
effective potential for the scalars and shows that the two conditions stated 
above are sufficient for the attractor phenomenon to hold. 

Second, we carry out numerical analysis. This requires a specific form of 
the effective potential, but allows us to go beyond the perturbative regime. 
The numerical analysis corroborates the perturbation theory results men- 
tioned above. In simple cases we have explored so far, we have found ev- 
idence for a only single basin of attraction, although multiple basins must 
exist in general as is already known from the SUSY cases. 

Finally, in some special cases, we solve the equations of motion exactly 
by mapping them a solvable Toda system. This allows us to study the black 
hole solutions in these special cases in some depth. Once again, in all the 
cases we have studied, we can establish the attractor phenomenon. 

It is straightforward to generalise these results to other settings. We 
find that the attractor phenomenon continues to hold in Anti-de Sitter space 
(AdS) and also in higher dimensions, as long as the two conditions men- 
tioned above are valid for a suitable defined effective potential. There is also 
possibly an attractor mechanism in de Sitter space (dS), but in the simplest 
of situations analysed here some additional caveats have to be introduced to 
deal with infrared divergences in the far past (or future) of dS space. 

This paper is structured as follows. Black holes in asymptotically flat 
four dimensional space are analysed first, in sections 2,3,4. The discussion is 



4 



extended to asymptotically flat space-times of higher dimension in section 5. 
Asymptotically AdS space is discussed next in section 6. 

As was mentioned above our analysis in the asymptotically flat and AdS 
cases is based on theories which have no potential for the scalars so that their 
values can vary at infinity. Some comments on this are contained in Section 
7. With Af ^ 1 SUSY such theories can arise, with the required couplings 
between scalars and gauge fields, and are at least technically natural. In the 
absence of supersymmetry there is no natural way to arrange this and our 
study is more in the nature of a mathematical investigation. We follow in 
section 8, with some comments on the attractor phenomenon in dS. Finally, 
in section 9 we show that non-extremal black holes do not have an attractor 
mechanism. Thus, the double-zero nature of the horizon is essential to draw 
the moduli to fixed values. 

Several important intermediate steps in the analysis are discussed in ap- 
pendices IATID1 

Some important questions are left for the future. First, we have not 
analysed the stability of these black hole solutions. It is unlikely that there 
are any instabilities at least in the S-wave sector. We do not attempt a 
general analysis of small fluctuations here. Second, in this paper we have 
not analysed string theory situations where such non-supersymmetric black 
holes can arise [Ej. This could include both critical and non-critical string 
theory. In case of M = 1 supersymmetry it would be interesting to explore 
if there is partial restoration of supersymmetry at the horizon. Given the 
rotational invariance of the solutions one can see that no supersymmetry is 
preserved in-between asymptotic infinity and the horizon in this case. 

Let us also briefly comment on some of the literature of especial relevance. 
The importance of the effective potential, V e ff, for M = 2 black holes was 
emphasised in jH]. Some comments pertaining to the non-supersymmetric 
case can be found for example in A similar analysis using an effec- 
tive one dimensional theory, and the Gauss-Bonett term, was carried out 
in [THl- Finally, while the thrust of the analysis is different, our results are 
quite closely related to those in [TH] which appeared while this paper was 
in preparation (see also (201 f° r the 3-dimensional case). In ^J] the entropy 
(including higher derivative corrections) is obtained from the gauge field La- 
grangian after carrying out a Legendre transformation with respect to the 
electric parameters. This is similar to our result which is based on V e ff. As 
was mentioned above, V e ff, is proportional to the electro-magnetic energy 
density i.e., the Hamiltonian density of the electro-magnetic fields, and is de- 
rived from the Lagrangian by doing a canonical transformation with respect 
to the gauge fields. For an action with only two-derivative terms, our results 



5 



and those in [TH] agree [2*Tj . 



2 Attractor in Four- Dimensional Asymptoti- 
cally Flat Space 

2.1 Equations of Motion 

In this section we consider gravity in four dimensions with U(l) gauge fields 
and scalars. The scalars are coupled to gauge fields with dilaton-like cou- 
plings. It is important for the discussion below that the scalars do not have 
a potential so that there is a moduli space obtained by varying their values. 

The action we start with has the form, 

S = ^J d'xV^GiR - 2(d<p i ) 2 - f ab (&)F« u F b n (1) 

Here the index i denotes the different scalars and a, b the different gauge fields 
and F® u stands for the field strength of the gauge field. fab(4>i) determines 
the gauge couplings, we can take it to be symmetric in a, b without loss of 
generality. 

The Lagrangian is 

C = (R- 2(d<p i ) 2 - f ah {4>i)F; v F h n (2) 

Varying the metric gives 3 , 

R»u - Id^idAi = 2f ab {<p i )F\ x F b u x + hs^C (3) 

The trace of the above equation implies 

R - 2(d(j) i ) 2 = (4) 

The equations of motion corresponding to the metric, dilaton and the gauge 
fields are then given by, 

R, v -2d^ i d v <j> i = f ab (cf> i )(2F a ^F b u x -lG flv F a KX F bKX ) (5) 

-±=d,(V^Gd^) = ~$(/ a6 )F%,F 6 »-' (6) 

d.jV^GfaMF'n = 0. 

3 In our notation G^ u refers to the components of the metric. 
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The Bianchi identity for the gauge field is, 

d„F vp + d u F Pfl + dpF^ = 0. (7) 

We now assume all quantities to be function of r. To begin, let us also 
consider the case where the gauge fields have only magnetic charge, general- 
isations to both electrically and magnetically charged cases will be discussed 
shortly. The metric and gauge fields can then be written as, 

ds 2 = -a{rfdt 2 + a(r)- 2 dr 2 + 6(r) W (8) 
F a = Q a m sin9d9 A# (9) 

Using the equations of motion we then get, 

Rtt = tfVefffa) (10) 

R9t = ^Veffifa) (11) 

where, 

Vefffc) = fa b (&)Q a m Q b m - (12) 

This function, V e ff, will play an important role in the subsequent discussion. 
We see from eq.(JTUJ) that up to an overall factor it is the energy density in 
the electromagnetic field. Note that V e ff(<pi) is actually a function of both 
the scalars and the charges carried by the black hole. 

2 

The relation, R tt = 75R00, after substituting the metric ansatz implies 
that, 

(a 2 (r)b 2 (r))" = 2. (13) 
The R rr — ^Rtt component of the Einstein equation gives 

j = -(dAf. (14) 
Also the Rrr component itself yields a first order "energy" constraint, 



- 1 + a 2 b' 2 + = -±(V eff m + a 2 b 2 ^) 2 (15) 

Finally, the equation of motion for the scalar 0j takes the form, 

d r {a 2 b 2 d r ^) = (16) 

We see that V e ff(<pi) plays the role of an "effective potential " for the scalar 
fields. 
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Let us now comment on the case of both electric and magnetic charges. 
In this case one should also include "axion" type couplings and the action 
takes the form, 

S = ^ J d*xV^G(R - 2(c^) 2 - Ui^F^F^ - \~f ah ^)F\ u F h pa e^°). 

(17) 

We note that f a b(4>i) is a function independent of fab(4>i), it can also be taken 
to be symmetric in a, b without loss of generality. 

The equation of motion for the metric which follows from this action is 
unchanged from eq. (p)|) . While the equations of motion for the dilaton and 
the gauge field now take the form, 

-A^e M (V=G^) = IW^F^F^" + Id^UF^F^t^ (18) 

d, (V=G (fMi)F h " V + \L b F b pa e^)) = 0. (19) 
With both electric and magnetic charges the gauge fields take the form, 
F a = f a \(t>i){Qe b ~ kQ c J^dt A dr + Q a m stn6d6 A #, (20) 

where Q a m , Q ea are constants that determine the magnetic and electric charges 
carried by the gauge field F a , and f ab is the inverse of f a j, 4 . It is easy to see 
that this solves the Bianchi identity eq.©, and the equation of motion for 
the gauge fields eq. (fT9|) . 

A little straightforward algebra shows that the Einstein equations for the 
metric and the equations of motion for the scalars take the same form as 
before, eq. lfEfl ITU Hoi Ho]) , with V e ff now being given by, 

Vefffa) = f a \Qea ~ facQ'JiQeb ~ ImQL) + fabQ^Ql- (21) 

As was already noted in the special case of only magnetic charges, V e ff is 
proportional to the energy density in the electromagnetic field and therefore 
has an immediate physical significance. It is invariant under duality trans- 
formations which transform the electric and magnetic fields to one-another. 

Our discussion below will use (fEfl ITU IT51 ITBjl and will apply to the general 
case of a black hole carrying both electric and magnetic charges. 

4 We assume that f a b is invertible. Since it is symmetric it is always diagonalisable. 
Zero eigenvalues correspond to gauge fields with vanishing kinetic energy terms, these can 
be omitted from the Lagrangian. 
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It is also worth mentioning that the equations of motion, eq. lfEfl fTH ITfij) 

above can be derived from a one-dimensional action, 

S=^Jdr Ua 2 b)'b' - aV(0') 2 - ■ (22) 

The constraint, eq. ([T5|) must be imposed in addition. 

One final comment before we proceed. The eq. (jl7|) can be further gener- 
alised to include non-trivial kinetic energy terms for the scalars of the form, 

J d^xV^G (- gij {fa)dfadfa>) . (23) 

The resulting equations are easily determined from the discussion above by 
now contracting the scalar derivative terms with the metric gij. The two 
conditions we obtain in the next section for the existence of an attractor are 
not altered due to these more general kinetic energy terms. 



2.2 Conditions for an Attractor 

We can now state the two conditions which are sufficient for the existence of 
an attractor. First, the charges should be such that the resulting effective 
potential, V e //, given by eq.(|2"T|). has a critical point. We denote the critical 
values for the scalars as fa = fa . So that, 

diV eff (fa ) = (24) 

Second, the matrix of second derivatives of the potential at the critical point, 

Mij = ^didjVefftta) (25) 

should have positive eigenvalues. Schematically we write, 

> (26) 

Once these two conditions hold, we show below that the attractor phe- 
nomenon results. The attractor values for the scalars are 5 fa = <piQ. 

The resulting horizon radius is given by, 

b 2 H = V eff (fao) (27) 

5 Scalars which do not enter in V e ff are not fixed by the requirement ea. (|24|l . The 
entropy of the extremal black hole is also independent of these scalars. 
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and the entropy is 



Sbh = \ A = vb 2 H . (2* 



There is one special solution which plays an important role in the discus- 
sion below. From eq.f|16j) we see that one can consistently set fa = <pio for all 
values of r. The resulting solution is an extremal Reissner Nordstrom (ERN) 
Black hole. It has a double-zero horizon. In this solution d r (j)i = 0, and a, b 
are 

ao(r) = (l - —), 6o(r) = r (29) 



where r# is the horizon radius. We see that a 2 ,, (a 2 ,)' vanish at the horizon 
while bo,b' are finite there. From eq. (jl5j) it follows then that the horizon 
radius bu is indeed given by 

rl = b 2 H = V eff (fa Q ), (30) 

and the black hole entropy is eq.(|55)l. 

If the scalar fields take values at asymptotic infinity which are small devi- 
ations from their attractor values we show below that a double-zero horizon 
black hole solution continues to exist. In this solution the scalars take the 
attractor values at the horizon, and a 2 , (a 2 )' vanish while b, b' continue to 
be finite there. From eq. (jl5j) it then follows that for this whole family of 
solutions the entropy is given by eq. ()28|) and in particular is independent of 
the asymptotic values of the scalars. 

For simple potentials V e ff we find only one critical point. In more com- 
plicated cases there can be multiple critical points which are attractors, each 
of these has a basin of attraction. 

One comment is worth making before moving on. A simple example of 
a system which exhibits the attractor behaviour consists of one scalar field 
coupled to two gauge fields with field strengths, F a ,a = 1,2. The scalar 
couples to the gauge fields with dilaton-like couplings, 

fM) = e aa4, 5 ah . (31) 

If only magnetic charges are turned on, 

V eff = e a ^(Q 1 ) 2 + e a ^(Q 2 ) 2 . (32) 

(We have suppressed the subscript "m" on the charges). For a critical point 
to exist ai and a 2 must have opposite sign. The resulting critical value of (j> 
is given by, 



= | | (33) 
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The second derivative, eq. ()25|) now is given by 



= -2a V (34) 



and is positive if a%, a>2 have opposite sign. 



This example will be useful for studying the behaviour of perturbation 
theory to higher orders and in the subsequent numerical analysis. 

As we will discuss further in section 7, a Lagrangian with dilaton-like 
couplings of the type in eq. (}3Tj) . and additional axionic terms ( which can be 
consistently set to zero if only magnetic charges are turned on), can always be 
embedded in a theory with M = 1 supersymmetry. But for generic values of 
a we do not expect to be able to embed it in an Af = 2 theory. The resulting 
extremal black hole, for generic a, will also then not be a BPS state. 



2.3 Comparison with the M = 2 Case 

It is useful to compare the discussion above with the special case of a BPS 
black hole in an M = 2 theory. The role of the effective potential, V e ff for 
this case was emphasised by Denef, j^j. It can be expressed in terms of a 
superpotential W and a Kahler potential K as follows: 

V eff = e K [K f W l W(D 1 W)* + \W\% (35) 

where DjW = djW + diKW. The attractor equations take the form, 

DiW = (36) 

And the resulting entropy is given by 

S BH = ir\W\ 2 e K . (37) 

with the superpotential evaluated at the attractor values. 

It is easy to see that if eq. (|36|) is met then the potential is also at a critical 
point, diV e ff = 0. A little more work also shows that all eigenvalues of the 
second derivative matrix, eq.(|25j) are also positive in this case. Thus the 
BPS attractor meets the two conditions mentioned above. We also note that 
from eq.(j35j) the value of V e ff at the attractor point is V e ff = e K \W\ 2 . The 
resulting black hole entropy eq. (|2*7l l2~%|) then agrees with eq.(pT7}). 

We now turn to a more detailed analysis of the attractor conditions below. 
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2.4 Perturbative Analysis 
2.4.1 A Summary 

The essential idea in the perturbative analysis is to start with the extremal 
RN black hole solution described above, obtained by setting the asymptotic 
values of the scalars equal to their critical values, and then examine what 
happens when the scalars take values at asymptotic infinity which are some- 
what different from their attractor values, fa = fao. 

We first study the scalar field equations to first order in the perturbation, 
in the ERN geometry without including backreaction. Let fa be a eigenmode 
of the second derivative matrix eq. (|23j) 6 . Then denoting, 5fa = fa — fa , 
neglecting the gravitational backreaction, and working to first order in 8 fa, 
we find that eq.(fTB|) takes the form, 

d r ((r-r H )%(8fa)) = S^, (38) 

where f3f is the relevant eigenvalue of ^didjV(fao). In the vicinity of the 
horizon, we can replace the factor 1/r 2 on the r.h.s by a constant and as we 
will see below, eq.([38|h has one solution that is well behaved and vanishes at 
the horizon provided f3f ^ 0. Asymptotically, as r — > oo, the effects of the 
gauge fields die away and eq. (J58j) reduces to that of a free field in flat space. 
This has two expected solutions, 8<pi ~ constant, and 8 fa ~ 1/r, both of 
which are well behaved. It is also easy to see that the second order differential 
equation is regular at all points in between the horizon and infinity. So once 
we choose the non-singular solution in the vicinity of the horizon it can be 
continued to infinity without blowing up. 

Next, we include the gravitational backreaction. The first order perturba- 
tions in the scalars source a second order change in the metric. The resulting 
equations for metric perturbations are regular between the horizon and infin- 
ity and the analysis near the horizon and at infinity shows that a double- zero 
horizon black hole solution continues to exist which is asymptotically flat 
after including the perturbations. 

6 More generally if the kinetic energy terms are more complicated, ea. (|23|l . these eigen- 
modes are obtained as follows. First, one uses the metric at the attractor point, gij(4>io), 
and calculates the kinetic energy terms. Then by diagonalising and rescaling one obtains 
a basis of canonically normalised scalars. The second derivatives of V e ff are calculated 
in this basis and gives rise to a symmetric matrix, ea. (|25[) . This is then diagonalised by 
an orthogonal transformation that keeps the kinetic energy terms in canonical form. The 
resulting eigenmodes are the ones of relevance here. 
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In short the two conditions, eq.(J2IJ), eq.(J21J), are enough to establish the 
attractor phenomenon to first non-trivial order in perturbation theory. 

In 4-dimensions, for an effective potential which can be expanded in a 
power series about its minimum, one can in principle solve for the pertur- 
bations analytically to all orders in perturbation theory. We illustrate this 
below for the simple case of dilaton-like couplings, eq.([31j). where the coeffi- 
cients that appear in the perturbation theory can be determined easily. One 
finds that the attractor mechanism works to all orders without conditions 
other than eq.fl2U), eq.flHJ) 7 . 

When we turn to other cases later in the paper, higher dimensional or 
AdS space etc., we will sometimes not have explicit solutions, but an analysis 
along the above lines in the near horizon and asymptotic regions and show- 
ing regularity in-between will suffice to show that a smoothly interpolating 
solution exists which connects the asymptotically flat region to the attractor 
geometry at horizon. 

To conclude, the key feature that leads to the attractor is the fact that 
both solutions to the linearised equation for 5(f) are well behaved as r — > oo, 
and one solution near the horizon is well behaved and vanishes. If one of these 
features fails the attractor mechanism typically does not work. For example, 
adding a mass term for the scalars results in one of the two solutions at 
infinity diverging. Now it is typically not possible to match the well behaved 
solution near the horizon to the well behaved one at infinity and this makes 
it impossible to turn on the dilaton perturbation in a non-singular fashion. 

We turn to a more detailed description of perturbation theory below. 
2.4.2 First Order Solution 

We start with first order perturbation theory. We can write, 

5(f)i = <j>i - 4> i0 = e(j)ii, (39) 

where e is the small parameter we use to organise the perturbation theory. 
The scalars <pi are chosen to be eigenvectors of the second derivative matrix, 
eq.flllj). 

From, eq.(JT3J), eq.(fri|). eq.(fT5j). we see that there are no first order cor- 
rections to the metric components, a, b. These receive a correction starting 

7 For some specific values of the exponent 7,-, ea. (|^T|l . though, we find that there can be 
an obstruction which prevents the solution from being extended to all orders. 
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at second order in e. The first order correction to the scalars (pi satisfies the 
equation, 

d r {albld r <f> iX ) = Ptfa. (40) 
% 

where, /3f is the eigenvalue for the matrix eq.f|25j) corresponding to the mode 
4>i. Substituting for a , b , from eq.(J2HJ) we find, 

f r _ r .„\i( ± v / i+^FA¥-i) 
0,i = c u [-^j (41) 

We are interested in a solution which does not blow up at the horizon, r = r H . 
This gives, 

<Pa = c u (^p L y\ (42) 

where 



* = !(\A^¥-i). («) 



Asymptotically, as r — > oo, 0a — > Cij, so the value of the scalars vary at 
infinity as is changed. However, since 7i > 0, we see from eq.p2|) that 4>n 
vanishes at the horizon and the value of the dilaton is fixed at 0jo regardless 
of its value at infinity. This shows that the attractor mechanism works to 
first order in perturbation theory. 

It is worth commenting that the attractor behaviour arises because the 
solution to eq. (}4T)j) which is non-singular at r = r#, also vanishes there. To 
examine this further we write eq. (|4U|) in standard form, [22j, 

^ + P(x)y + Q(x)y = 0, (44) 

with x — r — th, y — 4>a- The vanishing non-singular solution arises because 
eq.pOjl has a single and double pole respectively for P(x) and Q(x), as x — > 0. 
This results in (|44|) having a scaling symmetry as x — >• and the solution 
goes like x l1 near the horizon. The residues at these poles are such that the 
resulting indical equation has one solution with exponent 7$ > 0. In contrast, 
in a non-extremal black hole background, the horizon is still a regular singular 
point for the first order perturbation equation, but Q(x) has only a single 
pole. It turns out that the resulting non-singular solution can go to any 
constant value at the horizon and does not vanish in general. 
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2.4.3 Second Order Solution 



The first order perturbation of the dilaton sources a second order correction 
in the metric. We turn to calculating this correction next. 

Let us write, 

b = b + e 2 b 2 (45) 

2 2 2 

a = a + e a 2 
b 2 = b 2 Q + 2e 2 b 2 b , 

where b and a Q are the zeroth order extremal Reissner Nordstrom solution 
eq.JH. 

Equation (JT3J) gives, 

a 2 b 2 = (r - r H ) 2 + d x r + d 2 . (46) 

The two integration constants, d\, d 2 can be determined by imposing bound- 
ary conditions. We are interested in extremal black hole solutions with van- 
ishing surface gravity. These should have a horizon where b is finite and a 2 
has a "double- zero" , i.e., both a 2 and its derivative (a 2 )' vanish. By a gauge 
choice we can always take the horizon to be at r = Th- Both d\ and d 2 then 
vanish. Substituting eq. (j43|) in the equation (fT3j) we get to second order in e, 

2a 2 b b 2 + b\a 2 = 0. (47) 

Substituting for a , b then determines, a 2 in terms of b 2 , 

a2 = _ 2 (l-l£) 2 ^. (48) 
V r / r 



From eq.()14|) we find next that, 

h(r) = - £ ^|^r (Z^pL) ^ + Al r + A 2 r H (49) 



2(2 7i 

Ai,A 2 are two integration constants. The two terms proportional to these 
integration constant solve the equations of motion for b 2 in the absence of the 
0(e) 2 source terms from the dilaton. This shows that the freedom associated 
with varying these constants is a gauge degree of freedom. We will set A\ = 
A 2 = below. Then, b 2 is, 

wo-E^ir^r (50) 
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It is easy to check that this solves the constraint eq. (jl5j) as well. 

To summarise, the metric components to second order in e are given by 
eq. ()45|) with ao,b Q being the extremal Reissner Nordstrom solution and the 
second order corrections being given in eq. p%j) and eq. (p)Uj) . Asymptotically, 
as r — > oo, b 2 — > c x r, and, a 2 — > —2 x c, so the solution continues to 
be asymptotically flat to this order. Since 7^ > we see from eq. lpIHl IBTH) 
that the second order corrections are well defined at the horizon. In fact 
since b 2 goes to zero at the horizon, a 2 vanishes at the horizon even faster 
than a double- zero. Thus the second order solution continues to be a double- 
zero horizon black hole with vanishing surface gravity. Since b 2 vanishes the 
horizon area does not change to second order in perturbation theory and is 
therefore independent of the asymptotic value of the dilaton. 

The scalars also gets a correction to second order in e. This can be calcu- 
lated in a way similar to the above analysis. We will discuss this correction 
along with higher order corrections, in one simple example, in the next sub- 
section. 

Before proceeding let us calculate the mass of the black hole to second 
order in e. It is convenient to define a new coordinate, 

V = Kr) (51) 

Expressing a 2 in terms of y one can read off the mass from the coefficient of 
the 1/y term as y — > 00, as is discussed in more detail in Appendix 1X1 This 
gives, 

M = ^ + e 2 ^^ (52) 

i 

where r H is the horizon radius given by ()3U|). Since 7, is positive, eq.()43|). 
we see that as e increases, with fixed charge, the mass of the black hole 
increases. The minimum mass black hole is the extremal RN black hole 
solution, eq. (j2Uj) . obtained by setting the asymptotic values of the scalars 
equal to their critical values. 

2.4.4 An Ansatz to All Orders 

Going to higher orders in perturbation theory is in principle straightforward. 
For concreteness we discuss the simple example, eq.fETT]). below. We show 
in this example that the form of the metric and dilaton can be obtained 
to all orders in perturbation theory analytically. We have not analysed the 
coefficients and resulting convergence of the perturbation theory in great 
detail. In a subsequent section we will numerically analyse this example 
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and find that even the leading order in perturbation theory approximates 
the exact answer quite well for a wide range of charges. This discussion 
can be generalised to other more complicated cases in a straightforward way, 
although we will not do so here. 

Let us begin by noting that eq. (|13)) can be solved in general to give, 

a 2 b 2 = (r - r H f + d x r + d 2 (53) 

As in the discussion after eq (pU|) we set di = d 2 = 0, since we are interested 
in extremal black holes. This gives, 

a 2 b 2 = (r _ Th) 2^ ( 54) 

where r# is the horizon radius given by eq. (|30j) . This can be used to deter- 
mine a in terms of b. 

Next we expand b, (p and a 2 in a power series in e, 

oo 

b = fo + ^ e "fo n (55) 

n=l 
oo 

= 0o + ^e n 0n (56) 

71=1 
OO 

+ X> nfl ™ ( 57 ) 



~ _ „2 

n=l 



a = a 



where &o> a o are given by eq. (j29|) and O is given by eq 

The ansatz which works to all orders is that the n th order terms in the 
above two equations take the form, 

Mr) = c n (^p^j 7 (58) 
b n (r) = d n r ( r -^Ji] \ (59) 



and, 

\ n.7+2 

r — th 



r 

where 7 is given by eqs.(jHH) and in this case takes the value, 



(60) 



7 = ~ (VI - 2 ai a 2 - 1) . (61) 
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The discussion in the previous two subsections is in agreement with this 
ansatz. We found b\ = 0, and from eq.(|5U|) we see that 62 is of the form 
eq.()59|). Also, we found ai = and from eq. (j48j) a 2 is of the form eq.(|6U)). 
And from eq.(|42|) we see that <p\ is of form eq. (|58|) . We will now verify 
that this ansatz consistently solves the equations of motion to all orders in 
e. The important point is that with the ansatz eq.()581 each term in the 
equations of motion of order e n has a functional dependence (^—jr^) 2 ' yn - This 
allows the equations to be solved consistently and the coefficients c n , d n to 
be determined. 

Let us illustrate this by calculating eg. From eq. flHj) and eq. ljHlj) we see 
that the equation of motion for can be written in the form, 

2b(r) 2 d r ((r-r H ) 2 d r <j ) ) = e^+Qfa (62) 

To 0(e 2 ) this gives, 

27 

{2c 2 (e a ^Q 2 a 2 - 4r^ 7 (l + 2 7 )) + e^Qfcfc?) = (63) 

Notice that the term ( r ~ rg ) 27 has factored out. Solving eq.()63p for c 2 we 
now get, 

C2 = ( 64 ) 

More generally, as discussed in Appendix^l working to the required order 
in e we can recursively find, c n ,d n ,e n . 

One more comment is worth making here. We see from eq. (|5()jl that 62 
blows up when when 7 = 1/2. Similarly we can see from eq. (|A.17J) that b n 
blows up when 7 = - for b n . So for the values, 7 = -, where n is an integer, 
our perturbative solution does not work. 

Let us summarise. We see in the simple example studied here that a 
solution to all orders in perturbation theory can be found, b, <ft and a 2 are 
given by eq.lfH^jl. eq.(jHHJ) an d eq. (j6T?j) with coefficients that can be determined 
as discussed in Appendix |X] In the solution, a 2 vanishes at r# so it is the 
horizon of the black hole. Moreover a 2 has a double- zero at r#, so the solution 
is an extremal black hole with vanishing surface gravity. One can also see 
that b n goes linearly with r as r — ► 00 so the solution is asymptotically flat to 
all orders. It is also easy to see that the solution is non-singular for r ^ r#. 
Finally, from eq.(|58p we see that <p n = 0, for all n > 0, so all corrections 
to the dilaton vanish at the horizon. Thus the attractor mechanism works 
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to all orders in perturbation theory. Since all corrections to b also vanish at 
the horizon we see that the entropy is uncorrected in perturbation theory. 
This is in agreement with the general argument given after eq. (|28|) . Note 
that no additional conditions had to be imposed, beyond eq.p H which 
were already appeared in the lower order discussion, to ensure the attractor 
behaviour 8 . 



3 Numerical Results 

There are two purposes behind the numerical work we describe in this section. 
First, to check how well perturbation theory works. Second, to see if the 
attractor behaviour persists, even when e, eq.([39|). is order unity or bigger 
so that the deviations at asymptotic infinity from the attractor values are 
big. We will confine ourselves here to the simple example introduced near 
eq. ()31|) , which was also discussed in the higher orders analysis in the previous 
subsection. 

In the numerical analysis it is important to impose the boundary con- 
ditions carefully. As was discussed above, the scalar has an unstable mode 
near the horizon. Generic boundary conditions imposed at r — > oo will there- 
fore not be numerically stable and will lead to a divergence. To avoid this 
problem we start the numerical integration from a point fj near the horizon. 
We see from eq. ()58| I39* jl that sufficiently close to the horizon the leading or- 
der perturbative corrections 9 becomes a good approximation. We use these 
leading order corrections to impose the boundary conditions near the horizon 
and then numerically integrate the exact equations, eq. (jl3H4j) . to obtain the 
solution for larger values of the radial coordinate. 

The numerical integration is done using the Runge-Kutta method. We 
characterise the nearness to the horizon by the parameter 



8 In our discussion of exact solutions in section 4 we will be interested in the case, 
cti = —ot2- From eci. (|64l [A..17|l we see that the expressions for ci and d% become, 



It follows that in the perturbation series for <p and b only the C2 n +i(odd) terms and 
^2n(even) terms are non- vanishing respectively. 

9 We take the 0(e) correction in the dilaton, ea. (|42|l . and the 0(e 2 ) correction in b, a 2 , 
ea. (|151E5)l . This consistently meets the constraint ea. lfH)|) to 0(e 2 ). 



Sr 



(67) 



c 2 = 
d 3 = 



(65) 
(66) 
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Plot of (f> comparing numerical and 1 st order perturbation result (a[=-a2 = l.7) 








Figure 1: Comparison of numerical integration of <p with 1 st order perturba- 
tion result. The upper graph is a close up of the lower one near the horizon. 
The perturbation result is denoted by a dashed line. We chose «i, — a 2 = 1-7, 
Qi = 3, Q 2 = 3, 5r = 2.3 x 1CT 8 and C\ in the range [— |, |] 
• 0o = O. 
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Plot of <p comparing numerical and 1 st order perturbation result (ct\ =-0-2=3.1) 
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Figure 2: Comparison of numerical integration of <p with 1 st order perturba- 
tion result. The upper graph is a close up of the lower one near the horizon. 
The perturbation result is denoted by a dashed line. We chose «i, — a 2 = 3.1, 
Qi = 2, Q 2 = 3, 5r = 2.9 x 10~ 8 and c\ is in the range [— ~, |] 
. 0o = 0.13. 
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where is the point at which we start the integration. c\ refers to the 
asymptotic value for the scalar, eq.()42jl. 

In figs. (jll2j) we compare the numerical and 1 st order correction. The 
numerical and perturbation results are denoted by solid and dashed lines 
respectively. We see good agreement even for large r. As expected, as we 
increase the asymptotic value of <f), which was the small parameter in our 
perturbation series, the agreement decreases. 

Note also that the resulting solutions turn out to be singularity free and 
asymptotically flat for a wide range of initial conditions. In this simple exam- 
ple there is only one critical point, eq. ()33j) . This however does not guarantee 
that the attractor mechanism works. It could have been for example that as 
the asymptotic value of the scalar becomes significantly different from the 
attractor value no double-zero horizon black hole is allowed and instead one 
obtains a singularity. We have found no evidence for this. Instead, at least 
for the range of asymptotic values for the scalars we scanned in the numeri- 
cal work, we find that the attractor mechanism works with attractor value, 

eq.flHi. 

It will be interesting to analyse this more completely, extending this work 
to cases where the effective potential is more complicated and several critical 
points are allowed. This should lead to multiple basins of attraction as has 
already been discussed in the super symmetric context in e.g., P ITU]. 



4 Exact Solutions 

In certain cases the equation of motion can be solved exactly (221 • In this 
section, we shall look at some solvable cases and confirm that the extremal 
solutions display attractor behaviour. In particular, we shall work in 4 di- 
mensions with one scalar and two gauge fields, taking V e // to be given by 
eq.flm, 

K// = e"^(Q 1 ) 2 + e^(g 2 ) 2 . (68) 
We find that at the horizon the scalar field relaxes to the attractor value (|33|) 

aiQt 

which is the critical point of V e ff and independent of the asymptotic value, 
0oo. Furthermore, the horizon area is also independent of (j)^ and, as pre- 
dicted in section l2~2l it is proportional to the effective potential evaluated at 
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the attractor point. It is given by 



Area = 4irb 2 H = 47rT4 // (0 o ) (70) 
= 47r*(Qi) 2 ^(Q 2 ) 2 ^ (71) 

where 

(72) 



ax) 1 «i 



is a numerical factor. It is worth noting that when a,\ = —a 2 , one just has 

Urea = 2vr | Q 1 Q 2 1 (73) 



Interestingly, the solvable cases we know correspond to 7 = 1, 2, 3 where 
7 is given by (jIHJ). The known solutions for 7 = 1, 2 are discussed in i23 3 and 
references therein (although they fixed 0^ = 0). We found a solution for 
7 = 3 and it appears as though one can find exact solutions as long as 7 is 
a positive integer. Details of how these solutions are obtained can be found 
in the references and appendix IE1 

For the cases we consider, the extremal solutions can be written in the 
following form 

2 —-c - 

Q2Y s^-^ 



ei — * = m w a (74) 

b 2 = 4((Qi/i)- Q2 (Q2/ 2 r) ai!Q2 ( 75 ) 

a 2 = p 2 /b 2 (76) 

where p = r — rn and the fi are polynomials in p to some fractional power. 
In general the fi depend on 0^ but they have the property 

fi I Horizon 1- (77) 

Substituting (J77J) into ()74|75jl . one sees that that at the horizon the scalar 
field takes on the attractor value (J(j9|) and the horizon area is given by (J71|) . 

Notice that, when a = |a»|, ()74I75|) simplify to 



1^2 



,,,, |Q2| ( fl\ '" 



^ = (7g) 

b 2 = 2|Q 1 |iQ 2 |(/ 1 / 2 ) (79) 
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iog(Vfi 2 /ei 



(r) for various values of ip^ (extremal bh, a=2) 




Figure 3: Attractor behaviour for the case 7 = 1; a\, —o>2 = 2 
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0.2 - 




Figure 4: Attractor behaviour for the case 7 = 2; ati, —ot2 = 
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4.1 Explicit Form of the fa 



In this section we present the form of the functions fa mainly to show that, 
although they depend on 0^ in a non trivial way, they all satisfy (|77jl which 
ensures that the attractor mechanism works. It is convenient to define 

q2 = e ai<ho02 ( no summa tion) (80) 

which are the effective U(l) charges as seen by an asymptotic observer. For 
the simplest case, 7 = 1, we have 

/ l = l+(gr 1 | ai |(4 + af)-^p (81) 
Taking 7 = 2 and a\ = — a 2 — 2a/3 one finds 

fi = (l + (QiQ 2 )-*(Q? + <? 2 ¥p + l(QiQiQ2)-h 2 Y (82) 

Finally for 7 = 3 and ai\ = 4 , a 2 = —6 we have 

f 1= (l- Qa 2 p + Yla\p 2 - 6a p 3 ) s (83) 

f 2 = (1 - f a 2 p + 24a 2 p 2 - (48a 3 - 12a )p 3 + (48a^ - 24a a 2 ) p 4 ) « (84) 

where cto and a 2 are non-trivial functions of Qi. Further details are discussed 
in section El and appendix El The scalar field solutions for 7 = 1 and 2 are 
illustrated in figs. El and |U respectively. 

4.2 Supersymmetry and the Exact Solutions 

As mentioned above, the first two cases (7 = 1,2) have been extensively 
studied in the literature. 

The SUSY of the extremal a± = — a 2 = 2 solution is discussed in |24J. 
They show that it is supersymmetric in the context of M = 4 SUGRA. It 
saturates the BPS bound and preserves | of the supersymmetry - ie. it has 
M = 1 SUSY. There are BPS black-holes in this context which carry only 
one U(l) charge and preserve 5 of the supersymmetry. The non-extremal 
blackholes are of course non-BPS. 

On the other hand, the extremal a% = — a 2 = 2a/3 blackhole is non- 
BPS [213 • It ar i ses i n the context of dimensionally reduced 5D Kaluza-Klein 
gravity [26 and is embeddable in M = 2 SUGRA. There however are BPS 
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black-holes in this context which carry only one U(l) charge and once again 
preserve \ of the supersymmetry |2Tj . 

We have not investigated the supersymmetry of the 7 = 3 solution, we 
expect that it is not a BPS solution in a supersymmetric theory. 



5 General Higher Dimensional Analysis 
5.1 The Set-Up 

It is straightforward to generalise our results above to higher dimensions. We 
start with an action of the form, 

S = ^ J d d x^(R - 2{d<Pi) 2 - f ab {<j>i)F a F h ) (85) 

Here the field strengths, F a are (d — 2) forms which are magnetic dual to 
2-form fields. 

We will be interested in solution which preserve a SO(d — 2) rotation 
symmetry. Assuming all quantities to be function of r, and taking the charges 
to be purely magnetic, the ansatz for the metric and gauge fields is 10 

ds 2 = -a(r) 2 dt 2 + a(r)- 2 dr 2 + b(r) 2 d{l 2 d _ 2 (86) 
F a = Q a sm d ~ 3 6 sin d - 4 • • - dd A d<j> A • ■ • (87) 
F a = Q a sm d ~ 3 6 sin d - 4 • • • dB A d</> A • • • (88) 



The equation of motion for the scalars is 

d r {a 2 b d - 2 d^) = {d -^T eff . (89) 
Here V e ff, the effective potential for the scalars, is given by 

Veff = f ab {^)Q a Q b - (90) 

10 Black hole which carry both electric and magnetic charges do not have an SO(d — 2) 
symmetry for general d and we only consider the magnetically charged case here. The 
analogue of the two-form in 4 dimensions is the d/2 form in d dimensions. In this case 
one can turn on both electric and magnetic charges consistent with SO(d/2) symmetry. 
We leave a discussion of this case and the more general case of p-forms in d dimensions 
for the future. 
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From the (R rr — ^R a ) component of the Einstein equation we get, 

= J^gfM (91) 

The R rr component gives the constraint, 

_( d _ 2 ){(d - 3) - ab'(2a'b + {d- 3)ab')} = 2fa 2 a 2 b 2 - ^V eff (fa) 

(92) 

In the analysis below we will use eq.(|89|) to solve for the scalars and 
then eq.(|91|) to solve for b. The constraint eq.(|92j) will be used in solving 
for a along with one extra relation, R tt = (d — 3)^-Rgg, as is explained in 
appendixO These equations (aside from the constraint) can be derived from 
a one- dimensional action 



S = ±Jdr ((d-3)(d-2)b d - 4 (l + a 2 b' 2 ) + (d-2)b d - 3 (a 2 )'b 

-2a 2 b d - 2 (d r <P) 2 -^V eff ) 



(93) 



As the analysis below shows if the potential has a critical point at fa = fao 
and all the eigenvalues of the second derivative matrix dijV((j)io) are positive 
then the attractor mechanism works in higher dimensions as well. 



5.2 Zeroth and First Order Analysis 

Our starting point is the case where the scalars take asymptotic values equal 
to their critical value, fa = fao. In this case it is consistent to set the scalars 
to be a constant, independent of r. The extremal Reissner Nordstrom black 
hole in d dimensions is then a solution of the resulting equations. This takes 
the form, 

oo(r) = (l - jf^ b (r) = r (94) 

where Th is the horizon radius. From the eq. ()92)1 evaluated at r# we obtain 
the relation, 

rf- 3) = (d-A)\V eff (fa ) (95) 

Thus the area of the horizon and the entropy of the black hole are determined 
by the value of V e ff(fao), as in the four-dimensional case. 



27 



Now, let us set up the first order perturbation in the scalar fields, 

0i = 4>io + e<M (96) 
The first order correction satisfies, 

drialb^dr^) = J^fa (97) 

where, Pf is the eigenvalue of the second derivative matrix (,d ~^' dijV e ff((f)io) 
corresponding to the mode fa. This equation has two solutions. If f3f > 
one of these solutions blows up while the other is well defined and goes to 
zero at the horizon. This second solution is the one we will be interested in. 
It is given by, 

fai = ca(l - rfT 3 /^ 3 ) 7 ' (98) 

where 7 is given by 

7, = \ (-1 + y/l + A/3fr^/(d-3)A (99) 



5.2.1 Second order calculations (Effects of backreaction) 

The first order perturbation in the scalars gives rise to a second order cor- 
rection for the metric components, a, b. We write, 

b(r) = b (r) +e%(r) (100) 
a{rf = a (r) 2 + e 2 a 2 (r) (101) 
b{rf = b (r) 2 + 2e 2 6 2 (r)6 (r) (102) 

where do, &o are given in eq. (}94|) . 

From ()9ip one can solve for the second order perturbation b 2 ( r )- For 
simplicity we consider the case of a single scalar field, 0. The solution is 
given by double-integration form, 

=>b 2 (r) = d l r + d 2 - 1 



2(d - 3) (d- 4)7(27 - l)rf 
(-(d-4)F 1 - 2 7 , ^ m d " 3 ] 

+ (2 7 - 1) (^) d " 3 F [Q, 1 - 2 7 , 2^; (^) d " 3 ]) , (103) 
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where c[ = 2(d — 3) 2 c 2 7 2 r^ / \d — 2), a positive definite constant, and F is 
Gauss's Hypergeometric function. More generally, for several scalar fields, b 2 
is obtained by summing over the contributions from each scalar field. The 
integration constants d\,d%, in eq. (|ll)3|) . can be fixed by coordinate trans- 
formations and requiring a double-zero horizon solution. We will choose 
coordinate so that the horizon is at r = r#, then as we will see shortly the 
extremality condition requires both di,d2 to vanish. As r — > r# we have 
from eq. p()3jl that 



Since 7 > 0, we see that b 2 vanishes at the horizon and thus the area and the 
entropy are uncorrected to second order. At large r, b 2 (r) oc 0(r) + 0(1) + 
0(r 7 ~ 2d ) so asymptotic behaviour is consistent with asymptotic flatness of 
the solution. 

The analysis for a 2 is discussed in more detail in appendix O In the 
vicinity of the horizon one finds that there is one non-singular solution which 
goes like, 02(7") — > C(r — r^)^ 27+2 \ This solution smoothly extends to r — > 00 
and asymptotically, as r — > 00, goes to a constant which is consistent with 
asymptotic flatness. 

Thus we see that the backreaction of the metric is finite and well behaved. 
A double-zero horizon black hole continues to exist to second order in per- 
turbation theory. It is asymptotically flat. The scalars in this solution at the 
horizon take their attractor values irrespective of their values at infinity. 

Finally, the analysis in principle can be extended to higher orders. Unlike 
four dimensions though an explicit solution for the higher order perturbations 
is not possible and we will not present such an higher order analysis here. 

We end with Fig 5. which illustrates the attractor behaviour in asymp- 
totically flat 4 + 1 dimensional space. This figure has been obtained for the 
example, eq. (}3T| l32|) . The parameter 5r is defined in eq. (l67j) . 



Next we turn to the case of Anti-de Sitter space in four dimensions. Our 
analysis will be completely analogous to the discussion above for the four and 
higher dimensional case and so we can afford to be somewhat brief below. 




(104) 



6 Attractor in AdS^ 
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Numerical plot of 0(r) in 4+ 1 dimensions 
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Figure 5: Numerical plot of 0(r) with ol\ = —0:2 = 2 for the an extremal 
black hole in 4 + 1 dimensions displaying attractor behaviour. 

The action in 4-dim. has the form 

S = ^J d'xV^GiR -2k- 2(d<p i ) 2 - U^F" - lUi^e^F^) 

(105) 

where A = —3/L 2 is the cosmological constant. For simplicity we will discuss 
the case with only one scalar field here. The generalisation to many scalars is 
immediate and along the lines of the discussion for asymptotically flat four- 
dimensional case. Also we take the coefficient of the scalar kinetic energy 
term to be field independent. 

For spherically symmetric solutions the metric takes the form, eq.([5|). 
The field strengths are given by eq. (|20j) . This gives rise to a one dimensional 
action 

S =LJ dr U- (a 2 b 2 )" - 2a 2 bb" - 2a 2 6 2 (a r 0) 2 - 2^f- + ^) , (106) 

where V e ff is given by eq. (|21j) . The equations of motion, which can be derived 
either from eq. ([106)1 or directly from the action, eq. ()105|) are now given by, 

d r {a% 2 dA) = il07) 

j = -(9 r 0) 2 , (108) 
which are unchanged from the flat four- dimensional case, and, 

(a 2 (r)b 2 (r))" = 2(1 - 2A6 2 ), (109) 



30 



n 2 'h 2 ' — 1 1h 2 

- 1 + a 2 b' 2 + ^~ = -^(Ve ff m + aV(^) 2 + ^, (110) 
where the last equation is the first order constraint. 



6.1 Zeroth and First Order Analysis for V 

The zeroth order solution is obtained by taking the asymptotic values of the 
scalar field to be its critical values, O such that diV e ff(4>o) = 0. 

The resulting metric is now the extremal Reissner Nordstrom black hole 
in AdS space, j2H], given by, 

a (rf = (r - r " )2(L2 + ff + 2r " r + r2 > (111) 
b {r) = r (112) 



The horizon radius r# is given by evaluating the constraint eq. (|110J) at 
the horizon, 

(L 2 r 2 H + 2r|) 

(113) 



The first order perturbation for the scalar satisfies the equation 

b 2 



dr^lftdrfa) = ^01 (114) 



where, 

B 2 = 

2 



1 = ^?K//(0o). (H5) 



This is difficult to solve explicitly. 

In the vicinity of the horizon the two solutions are given by 

4>i = C±[r-T H ) t± (H6) 

If V"ff(<fio) > one of the two solutions vanishes at the horizon. We are 
interested in this solution. It corresponds to the choice, 

<h = C{T-r H y, (117) 

where, 

7= V »' , (H8) 
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and, 5 = - — L2 ' . As discussed in the appendix |D] this solution behaves at 
r — > oo as 0i — > C\ + C 2 jr i . Also, all other values of r, besides the horizon 
and oo, are ordinary points of the second order equation eq. ()114|) . All this 
establishes that there is one well-behaved solution for the first order scalar 
perturbation. In the vicinity of the horizon it takes the form eq. (|116|) with 
eq. ()118p . and vanishes at the horizon. It is non-singular everywhere between 
the horizon and infinity and it goes to a constant asymptotically at r — > oo. 

We consider metric corrections next. These arise at second order. We 
define the second order perturbations as in eq.(|45|). The equation for b 2 from 
the second order terms in eq. (jl08|) takes the form, 

6£ = -r(0Ur)) 2 , (119) 

and can be solved to give, 

b 2 (r) = - I I [r(0Ur)) 2 ] (120) 

Jr H Jr H 

We fix the integration constants by taking take the lower limit of both inte- 
grals to be the horizon. We will see that this choice gives rise to an double- 
zero horizon solution. Since is well behaved for all r# < r < 00 the 
integrand above is well behaved as well. Using eq. ()116j) we find that in the 
near horizon region 

b 2 ~ (r-r H ) {27) (121) 

At r — > 00 using the fact that 0i — ► C\ + C2/V 3 we find 

b 2 ~ D x r + D 2 + D 3 /r 6 . (122) 
This is consistent with an asymptotically AdS solution. 

Finally we turn to a 2 - As we show in appendix [Dl a solution can be found 
for a 2 with the following properties. In the vicinity of the horizon it goes 
like, 

a 2 oc(r-r H )^ +2 \ (123) 

and vanishes faster than a double- zero. As r — > 00, a 2 — > d\r and grows more 
slowly than a 2 ,. And for r# < r < 00 it is well-behaved and non-singular. 

This establishes that after including the backreaction of the metric we 
have a non-singular, double-zero horizon black hole which is asymptotically 
AdS. The scalar takes a fixed value at the horizon of the black hole and the 
entropy of the black hole is unchanged as the asymptotic value of the scalar 
is varied. 
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Let us end with two remarks. In the AdS case one can hope that there 
is a dual description for the attractor phenomenon. Since the asymptotic 
value of the scalar is changing we are turning on a operator in the dual 
theory with a varying value for the coupling constant. The fact that the 
entropy, for fixed charge, does not change means that the number of ground 
states in the resulting family of dual theories is the same. This would be 
worth understanding in the dual description better. Finally, we expect this 
analysis to generalise in a straightforward manner to the AdS space in higher 
dimensions as well. 

Fig. 6 illustrates the attractor mechanism in asymptotically AdS 4 space. 
This Figure is for the example, eq. (j31| E2J)- The cosmological constant is 
taken to be, A = —2.91723, in k — 1 units. 



Numerical plot of 0(r) in AdS4 





Gl=-7=r fi2=V2 

V2 



=0.346574 <Sr=0.0150842 A=-2.91723 



Figure 6: Numerical plot of <p(r) with a.\ = —a^ = 2 for the an extremal 
black hole in AdSi displaying attractor behaviour. 



7 Additional Comments 



The theories we considered in the discussion of asymptotically flat space- 
times and AdS spacetimes have no potential for the scalars. We comment on 
this further here. 

Let us consider a theory with M = 1 supersymmetry containing chiral 
superfields whose lowest component scalars are, 

S i = (f) i + ia i (124) 
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We take these scalars to be uncharged under the gauge symmetries. These 
can be coupled to the superfields by a coupling 



Such a coupling reproduces the gauge kinetic energy terms in and eq. (jl()5|) . 
eq. (|106|) . (we now include both 4>u °>i m the set of scalar fields which we 
denoted by 0j in the previous sections). 

An additional potential for the scalars would arise due to F-term con- 
tributions from a superpotential. If the superpotential is absent we get the 
required feature of no potential for these scalar. Setting the superpotential 
to be zero is at least technically natural due to its non-renormalisability. 

In a theory with no supersymmetry there is no natural way to suppress 
a potential for the scalars and it would arise due to quantum effects even 
if it is absent at tree-level. In this case we have no good argument for not 
including a potential for the scalar and our analysis is more in the nature of 
a mathematical investigation. 

The absence of a potential is important also for avoiding no-hair theo- 
rems which often forbid any scalar fields from being excited in black hole 
backgrounds j2H] • In the presence of a mass m in asymptotically flat four di- 
mensional space the two solutions for first order perturbation at asymptotic 
infinity go like, 



We see that one of the solutions blows up as r — > oo. Since one solution to the 
equation of motion also blows up in the vicinity of the horizon, as discussed 
in section 2, there will generically be no non-singular solution in first order 
perturbation theory. This argument is a simple-minded way of understanding 
the absence of scalar hair for extremal black holes under discussion here. In 
the absence of mass terms, as was discussed in section 2, the two solutions at 
asymptotic infinity go like ~ const and <f> ~ 1/r respectively and are both 
acceptable. This is why one can turn on scalar hair. The possibility of scalar 
hair for a massless scalar is of course well known. See ["JU|, [231, for some 
early examples of solutions with scalar hair, [3T| 13*2*1 "3*3*1 13*4*1 , for theorems on 
uniqueness in the presence of such hair, and [S] for a discussion of resulting 
thermodynamics. 

In asymptotic AdS space the analysis is different. Now the (mass) 2 for 
scalars can be negative as long as it is bigger than the BF bound. In this 
case both solutions at asymptotic infinity decay and are acceptable. Thus, 
as for the massless case, it should be possible to turn on scalar fields even 




(125) 



~ C ie m 7r, ~ C 2 e- mr /r. 



(126) 
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in the presence of these mass terms and study the resulting black holes solu- 
tions. Unfortunately, the resulting equations are quite intractable. For small 
(mass) 2 we expect the attractor mechanism to continue to work. 

If the (mass) 2 is positive one of the solutions in the asymptotic region 
blows up and the situation is analogous to the case of a massive scalar in flat 
space discussed above. In this case one could work with AdS space which 
is cut off at large r (in the infrared) and study the attractor phenomenon. 
Alternatively, after incorporating back reaction, one might get a non-singular 
geometry which departs from AdS in the IR and then analyse black holes 
in this resulting geometry. In the dual field theory a positive (mass) 2 cor- 
responds to an irrelevant operator. The growing mode in the bulk is the 
non-normalisable one and corresponds to turning on a operator in the dual 
theory which grows in the UV. Cutting off AdS space means working with 
a cut-off effective theory. Incorporating the back-reaction means finding a 
UV completion of the cut-off theory. And the attractor mechanism means 
that the number of ground states at fixed charge is the same regardless of 
the value of the coupling constant for this operator. 



8 Asymptotic de Sitter Space 

In de Sitter space the simplest way to obtain a double-zero horizon is to 
take a Schwarzschild black hole and adjust the mass so that the de Sitter 
horizon and the Schwarzschild horizon coincide. The resulting black hole 
is the extreme Schwarzschild-de Sitter spacetime jHHj- We will analyse the 
attractor behaviour of this black hole below. The analysis simplifies in 5- 
dimensions and we will consider that similar analysis can be carried 

out in other dimensions as well. Since no charges are needed we set all the 
gauge fields to zero and work only with a theory of gravity and scalars. Of 
course by turning on gauge charges one can get other double-zero horizon 
black holes in dS, their analysis is left for the future. 

We start with the action of the form, 

S=-^ J d 5 xy/^G(R - 2(d<P) 2 - V(<P)) (127) 

Notice that the action now includes a potential for the scalar, V(<fi), it will 
play the role of V e ff in our discussion of asymptotic flat space and AdS space. 
The required conditions for an attractor in the dS case will be stated in terms 
of V. A concrete example of a potential meeting the required conditions will 
be given at the end of the section. For simplicity we have taken only one 
scalar, the analysis is easily extended for additional scalars. 
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The first condition on V is that it has a critical point, V'(<po) = 0. We 
will also require that V((j) ) > 0. Now if the asymptotic value of the scalar 
is equal to its critical value, <fio, we can consistently set it to this value for 
all times t. The resulting equations have a extremal black hole solution 
mentioned above. This takes the form 

ds 2 = - {t2/L t2 _ L/2)2 df + (tVL ; 2 L/2) V + tHnl (128) 

Notice that it is explicitly time dependent. L is a length related to V(4>o) 
by , V((f>o) = |§. And t = ±-^= is the location of the double-zero horizon. 
A suitable near-horizon limit of this geometry is called the Nariai solution, 



8.1 Perturbation Theory 

Starting from this solution we vary the asymptotic value of the scalar. We 
take the boundary at t — > — oo as the initial data slice and investigate what 
happens when the scalar takes a value different from 0o a s t — > — oo. Our 
discussion will involve part of the space-time, covered by the coordinates 
in eq. (|128|) . with — oo < t < tn = — We carry out the analysis in 
perturbation theory below. 

Define the first order perturbation for the scalar by, 

= 0o + e0i 

This satisfies the equation, 

dtialbldtfa) = jK"(0 o )0i (129) 

where a = - ^ L ~ L ^ 2 \ b = t. This equation is difficult to solve in general. 
In the vicinity of the horizon t = tn, we have two solutions which go like, 

<!> 1 = C±(t-t H ) =2± ^ (130) 

where 

k 2 = (131) 

We see that one of the two solutions in eq. ()13U|) is non-divergent and in fact 
vanishes at the horizon if 

V"(</> ) < 0. (132) 
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We will henceforth assume that the potential meets this condition. Notice 
this condition has a sign opposite to what was obtained for the asymptotically 
flat or AdS cases. This reversal of sign is due to the exchange of space and 
time in the dS case. 

In the vicinity of t — ► — oo there are two solutions to eq. (|129j) which go 
like, 

fa = C ± \t\ p± (133) 

where 

p± = 2(-l ± V1 + KV4). (134) 

If the potential meets the condition, eq. f|132|) then k 2 > and we see that 
one of the modes blows up at t — > — oo. 



8.2 Some Speculative Remarks 

In view of the diverging mode at large \t\ one needs to work with a cutoff 
version of dS space 11 . With such a cutoff at large negative t we see that 
there is a one parameter family of solutions in which the scalar takes a fixed 
value at the horizon. The one parameter family is obtained by starting with 
the appropriate linear combination of the two solutions at t — > — oo which 
match to the well behaved solution in the vicinity of the horizon. While we 
will not discuss the metric perturbations and scalar perturbations at second 
order these too have a non-singular solution which preserves the double-zero 
nature of the horizon. The metric perturbations also grow at the boundary 
in response to the growing scalar mode and again the cut-off is necessary to 
regulate this growth. This suggests that in the cut-off version of dS space 
one has an attractor phenomenon. Whether such a cut-off makes physical 
sense and can be implemented appropriately are question we will not explore 
further here. 

One intriguing possibility is that quantum effects implement such a cut- 
off and cure the infra-red divergence. The condition on the potential eq. ()132|) 
means that the scalar has a negative (mass) 2 and is tachyonic. In dS space 
we know that a tachyonic scalar can have its behaviour drastically altered 
due to quantum effects if it has a (mass) 2 < H 2 where H is the Hubble 
scale of dS space. This can certainly be arranged consistent with the other 
conditions on the potential as we will see below. In this case the tachyon can 
be prevented from "falling down" at large \t\ due to quantum effects and the 
infrared divergences can be arrested by the finite temperature fluctuations of 

11 This is related to some comments made in the previous section in the positive (mass) 2 
case in AdS space. 
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dS space. It is unclear though if any version of of the attractor phenomenon 
survives once these quantum effects became important. 

We end by discussing one example of a potential which meets the various 
conditions imposed above. Consider a potential for the scalar, 

V = A ie Ql * + A 2 e a2<t> . (135) 

We require that it has a critical point at <fi = 0o and that the value of the 
potential at the critical point is positive. The critical point for the potential 
eq. (|135|) is at, 

- 1 - (136) 



Requiring that V((f) ) > tells us that 



V^o) = A 2 e a ^° - > (137) 

Finally we need that V"{<po) < this leads to the condition, 

V"(<f>o) = A 2 e a ^°a 2 (a 2 - a t ) < (138) 

These conditions can all be met by taking both a±, a 2 > 0, a 2 < a±, A 2 > 
and Ai < 0. In addition if a 2 a\ 3> 1 the resulting —{mass) 2 3> H 2 . 



9 Non-Extremal = Unattractive 

We end the paper by examining the case of an non-extremal black hole which 
has a single-zero horizon. As we will see there is no attractor mechanism in 
this case. Thus the existence of a double-zero horizon is crucial for the 
attractor mechanism to work. 

Our starting point is the four dimensional theory considered in section 
2 with action eq.(|17j). For simplicity we consider only one scalar field. We 
again start by consistently setting this scalar equal to its critical value, <po, 
for all values of r, but now do not consider the extremal Reissner Nordstrom 
black hole. Instead we consider the non-extremal black hole which also solves 
the resulting equations. This is given by a metric of the form, eq.(JSJ), with 

a 2 (r) = (l - r -±) (l - r f) , b(r) = r (139) 

where r± are not equal. We take r + > r_ so that r + is the outer horizon 
which will be of interest to us. 
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The first order perturbation of the scalar field satisfies the equation, 

<9 r (aVt^) = ^^01 (140) 

In the vicinity of the horizon r = r + this takes the form, 

dyiydyfa) = afa (141) 

where a is a constant dependent on V"((f>o),r + , r_, and y = r — r+. 

This equation has one non-singular solution which goes like, 

0i = C + C iy + --- (142) 

where the ellipses indicate higher terms in the power series expansion of <f>i 
around y — 0. The coefficients Ci, C2, ■ • • are all determined in terms of Co 
which can take any value. Thus we see that unlike the case of the double- 
horizon extremal black hole, here the solution which is well-behaved in the 
vicinity of the horizon does not vanish. 



Asymptotically, as r — > 00 both solutions to eq. ()140|) are well defined and 
go like 1/r, constant respectively. It is then straightforward to see that one 
can choose an appropriate linear combination of the two solutions at infinity 
and match to the solution, eq. ()142j) in the vicinity of the horizon. The impor- 
tant difference here is that the value of the constant Co in eq. ()142j) depends 
on the asymptotic values of the scalar at infinity and therefore the value 
of 4> does not go to a fixed value at the horizon. The metric perturbations 
sourced by the scalar perturbation can also be analysed and are non-singular. 
In summary, we find a family of non-singular black hole solutions for which 
the scalar field takes varying values at infinity. The crucial difference is that 
here the scalar takes a value at the horizon which depends on its value at 
asymptotic infinity. The entropy and mass for these solutions also depends 
on the asymptotic value of the scalar 12 . 

It is also worth examining this issue in a non-extremal black holes for an 
exactly solvable case. 

If we consider the case |aj| = 2, sectionHJ the non-extremal solution takes 
on a relatively simple form. It can be written [2*4*] 

exp(20) = e 2 ^^±H 
(r-S) 



12 An intuitive argument was given in the introduction in support of the attractor mech- 
anism. Namely, that the degeneracy of states cannot vary continuously. This argument 
only applies to the ground states. A non-extremal black hole corresponds to excited states. 
Changing the asymptotic values of the scalars also changes the total mass and hence the 
entropy in this case. 
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where 13 

r±=M±r r = ^ M 2 + E 2 - Q\ - Q\. (144) 
and the Hamiltonian constraint becomes 

Y? + M 2 -Ql-Ql = h y r + -r_y. (145) 

The scalar charge, E, defined by (f> ~ + — , is not an independent param- 
eter. It is given by 

E = (M6) 

There are horizons at r = r±, the curvature singularity occurs at r = E and 
ro characterises the deviation from extremality. We see that the non-extremal 
solution does not display attractor behaviour. 

Fig. [7| shows the behaviour of the scalar field 14 as we vary (p^ keeping 
M and Qi fixed. The location of the horizon as a function of r depends on 
0oo, eq. (|144|) . The horizon as a function of 0^ is denoted by the dotted line. 
The plot is terminated at the horizon. 

In contrast, for the extremal black hole, 

M= m±m bJ&i-i&i, (147) 

so ()143j) gives 

M-E IQil' U48j 

which is indeed the attractor value. 
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Figure 7: Plot 0(r) with ot\ = — a 2 = 2 for the non-extremal black hole with 
M, Qi held fixed while varying 0^ . The dotted line denotes the outer horizon 
at which we terminate the plot. It is clearly unattractive. 

A Perturbation Analysis 



A.l Mass 

Here, we first calculate the mass of the extremal black hole discussed in 
section 2.2. From eq.([5UJ). for large r, 6 2 is given by, 

b 2 = C r + d (A.l) 

where 

° h (A.2) 



2(2 7 -l) 

' - 0h (A ' 3) 

Now, we can easily write down the expression for a 2 using eq. (|48jl . We choose 
coordinate y as introduced in eq.(jHH) such that at large r, 

r 2 + 2e 2 (cr 2 + dr) = y 2 (A.4) 

- = -(l + 6 2 (c + ^)) (A.5) 

r y y 

We use the extremality condition (J53|) to find, 

a(r) = (A.6) 
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Using, eq. ljA.ll IA.6|) one finds that asymptotically, as r —>■ oo the metric 
takes the form, 

2{r H + e 2 (cr H + d))\ „ 1 2 2 2 



rfs = " ^ " 1 " ) dt + (1 ^gB ) ^ + f 

(A.7) 

where t is obtained by rescaling t and <if2 2 denotes the metric of S 2 . The 
mass M of the black hole is then given by the 1/y term in the g yy component 
of the metric. This gives, 

M = r H + e 2r -^ 1 . (A.8) 



A. 2 Perturbation Series to All Orders 

Next we go on to discuss the perturbation series to all orders, Using (|55|) for 
b and (JHHJ) for (f> in eq.((T3J)and eq.([H2|). we get, 

^-EE^^-h ( A -9) 

i=0 j=0 



26 i 6 fc _ i _ i ((r - r^) 2 ^)' = Q^o^* (A.10) 

where 

Vjfc = V" gi J 2 / i ar i+n2+ - + " fc . (A.11) 
' n x \n 2 \ ...n k \ 

{"l.»2 — n fc> 
J] mn m =k 

(A.12) 

After substituting our ansatz (Jo"%|) and (|59p. the above equations give, 

^(A/y - 1)4 = -7 ^ 3{k-i- j)diCjC k , (A. 13) 

i+j<k 

and 

fc(fc7 + l)c fc + T fc = ( 1 + l)( Ck + S k ) (A.14) 
where and T k are given by 

"1 c 2 • • • c fc-l 

ni!n 2 ! . ..n k -i ] . V" 1 



S fc = £ Cl , C2 ,--- Cfc - 1 , (aF"'- 1 + aF"'- 1 ) (A.15) 



{n 1 ,n 2 ...n k _ 1 } 
Y, mn m =k 
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and 

Tk=J2 l ^ + WA-i-fr- (A.16) 



Kk 



Then solving for dk and Ck gives 



1 \ ~^ 

dk = ~ k(k -1) ^ ^ k ~ i ~ j) d i e 3 e k-i-3 ( A - 17 ) 
^ ^ ' i+j<k 

(7 + l)gfc - T k 

k ~ ((fc + l) 7 + l)(ifc-l) [ } 

Finally, can be obtained using eq. (|54|) . eq.(|59j). It can be verified that the 
ansatz, eq. (|551 EH EOJ) with the coefficients eq. (|A.l 71 IA.I8|) also solves the 
constraint eq.(fT5|). 



B Exact Analysis 

Exact solutions can be found by writing the equations of motion as gener- 
alised Toda equations |3*7] . which may, in certain special cases, be solved 
exactly [23J - we rederive this result in slightly different notation below. As 
noted in in a marginally different context, the extremal solutions, are, in 
appropriate variables, polynomial solutions of the Toda equations. The poly- 
nomial solutions are much easier to find and are related to the functions /j 
mentioned in section For ease of comparison we occasionally use notation 
similar to 



B.l New Variables 

To recast the equations of motion into a generalised Toda equation we define 
the following new variables 

ui = (ft u 2 = log a z = logafe • = d T — a 2 b 2 d r (B.l) 

In terms of r, r is given by 

dr 1 ( r — r+ \ /rx s 

log ± (B.2) 



a 2 b 2 (r + — r_) \r — T- r 
where r± are the integration constants of ()13p. In general (JT3~j) implies 

a 2 b 2 = (r -r+)(r -r_). (B.3) 
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Notice that 

t — ► as r — > oo (B.4) 
r — > — oo as r — * r + (B.5) 

When we have a double- zero horizon, r# = r±, r takes the simple form 

r -1 = -(r-r-tf). (B.6) 

Since we are mainly interested in solutions with double- zero horizons, in what 
follows it will be convenient to work with a new radial coordinate, p, defined 
by 

p = -r-\ (B.7) 
which has the convenient property that pu = 0. 



B.2 Equivalent To da System 

In terms of these new variables the equations of motion become 

ut = \a x e 2u2 + a ^Ql + \oc 2 e 2u * +a ^Ql (B.8) 
u 2 = e 2u2+aiUl Ql + e 2u2+a2Ul Q 2 2 (B.9) 
z = e 2z (B.10) 

u 2 + u 2 -z 2 + e 2z - e 2u2+aiUl Ql - e 2u2+a2Ul Q 2 2 = (B.ll) 

(IB.10J) decouples from the other equations and is equivalent to (fT3|) . Finally 
making the coordinate change 

Xi = n- l Uj + trig log ((«i - a 2 )Q 2 ) (B.12) 

where 

■n- 1 = ( [ 

, -2 ai 
and 

= ^77 V ( 4 + ( B - 14 ) 

we obtain the generalised 2 body Toda equation 

.Y, f "'" v . (B.15) 

together with 

^ (f V,/,/ (; ,V ; - e^i) = (ai - a 2 )£ (B.16) 
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= ( "" 2 ) (B.13) 



where £ = \{r + 
given by 



— r_) 2 . After solving the above, the original fields will be 



e 



(ai— a.2) 



Q2 e l {aiXl+ a 2 X 2 ) (B1?) 



Ql ^ 

a 2 = e ^ {Xl+X ' 2) /<> (B.18) 
b 2 BP ( r _ r+ )( r _ r _)/ a 2 (B.19) 

where 



0=(ai-a2)Qr i -° a Q 2 ai - aa ( B - 20 ) 

B.3 Solutions 

B.3.1 Case I: 7 = 1 <^> a x a 2 = -4 

In this case, is diagonal 

m = diag(«i/2, -a 2 /2), (B.21) 
so the equations of motion decouple: 

Xi = e' Xi . (B.22) 



()B.22|) has solutions 



2 , ( 4c * 
,— r log 1 i . -1 2 / / 
|CKj| \ I «j I smh (Cj(r — dj J J 



^ = TTT lo S ( — T27V JT^ J (B.23) 



The integration constants are fixed by imposing asymptotic boundary con- 
ditions and requiring that the solution is finite at the horizon. Letting 

Fi = sinh( Ci (r - di)) (B.24) 

in terms of <fi and a we get 

3 (Ql-Q 2 )0 



^2 p |(aiXi+a2X 2 ) 
"2 A / Q2F2C1 



Oil J \Q1F1C2 

a 2 = e ^ (M2) /0 



(B.25) 



o ~ Q 2 2 - 



QiFj \Q 2 F 
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As r — > r + (ie. r — > — oo) the scalar field goes like 

e (ai-a a )* _ e 2( Cl -c 2 )r ^g^gj 

so we require 

c := ci = c 2 (B.27) 
for a finite solution at the horizon. Also at the horizon 

b 2 ~ (r - r+)/a 2 ~ e ((»-+-r-)-2c)r (R28 ) 

which necessitates 

(r+ - r_) = 2c (B.29) 
To find the extremal solutions we take the limit c — > which gives 

(a!-*)* = (_0»\W2f (B3()) 

& 2 = 4(Qi/i)^(Q 2 p 2 )^ (B.31) 

a 2 = p 2 /6 2 (B.32) 

where 

fi = 1 + d<p. (B.33) 
Requiring — > 0^ and a — > 1 as r — ► oo fixes 



di = Qi \ (B.34) 
y «i — a 2 

where as before 

Q2 = e *<fiooQ2 ( no summa tion). (B.35) 
For comparison with the non-extremal solution in this case see section |0 



B.3.2 Case II: 7 = 2 and oc x = -a 2 = 2^ 

In this case, becomes 

m=( ~J\ ■ ( B - 36 ) 



V3 a/3 / 



It is convenient to use the coordinates 



& = -=Xi - V3 log v 7 ^ (B.37) 
v3 
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so the equations of motion are the two particle Toda equations 

q \ = e 2qi - q2 (B.38) 
q 2 = e 2q2 - qi . (B.39) 

These maybe integrated exactly but the explicit form is, in general, a little 
complicated. Fortunately we are mainly interested in extremal solutions 
which have a simpler form [3Hj. As in, [HH], taking the ansatz that e~ qi is a 
second order polynomial one finds 

e~ qi = ao + ^r + K 2 (B.40) 

e"" 2 = a \ - a + a X T + 1 -t 2 (B.41) 

Finally, returning to the original variables and imposing the asymptotic 
boundary conditions gives the solution 

^ - (t)W (b - 42) 

b 2 = 2Q 1 Q 2 f 1 f 2 (B.43) 
a 2 = p 2 /b 2 (BAA) 



where 

fi= (i + {QiQ2)-HQl + Q2hh+UQiQiQi)~~ 3 P 2 y ( R45 ) 

as quoted in section HI 

For completeness we note that the general, non-extremal solution of [2E1 
l3*U] . modified for a non-zero asymptotic value of 4>, is 

exp(40/v / 3) =e^/^ — (B.46) 

Pi 

a 2 = (r ~ r+)(r - r - } (B.47) 

VP1P2 

b 2 = (B.48) 



where 



Pi = (r - r i+ )(r - nJ) (B.49) 



2 - AY, 

n± = t rS ± QiJ — — (B.50) 

and scalar charge, S, which is again not an independent parameter, is given 

by - 2 - 2 

7! = 2M(A-1) + 2M(A + 1) TIM ( ' 51) 
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B.3.3 Case III: 7 = 3 and on = 4 a 2 = -6 



In this case, m^- becomes 



m= ( _\ 2 X } (B.52) 



Making the coordinate change 



qi = |Xi-bg2 (B.53) 
g 2 = X 2 -log2 (B.54) 



The equations of motion are 



q\ = e 291 "" 2 (B.55) 
g 2 = e 292 - 292 (B.56) 



Now consider the three particle Toda system 



q\ = e 2qi ' q2 (B.57) 
q 2 = e ^2-<n-i3 (B.58) 
q 3 = e 2 « (B.59) 

which may be integrated exactly. Notice that by identifying q± and qs we 
obtain (jB.57HB.59j) . Once again the general solution is slightly complicated 
but taking the ansatz that e~ qi is a polynomial one finds 

e -« = a + 1a\r + a 2 r 2 + -r 3 (B.60) 

6 

2 1 

e~ q2 = Aa\ - 2a a 2 + (4a 3 , - a )r + 2a 2 r 2 + -a 2 r 3 + — r 4 (B.61) 



3 12 



Rewriting in terms of the original fields we get 



100 _ / 



2 



) exp (2Xi - 3X 2 ) (B.62) 



6 (()>Y fh x 12 



4 VQi/ \/i 



(B.63) 



p 2 10gf Qj exp f- jU a - ix 2 ) (B.64) 
5 /2\* * 



, } . 8 , QfQI/1/2 (B.65) 
a 2 = p 2 /6 2 (B.66) 
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where 

/i = (1 - 6a 2 p + \2a\p 2 - 6a p 3 ) 1 (B.67) 

h = (1 - f a 2 p + 24a 2 p 2 - (48^ - 12a )p 3 + (48a* - 24a a 2 ) p 4 )^ 

(B.68) 

At the horizon we do indeed have <fi at the critical point of V e ff\ 



~2Qi 



and b 2 given by \/ e/ /(0 o ): 



2 V3; " " VQi 

Imposing the asymptotic boundary conditions we get 



2f „ _ 2^ 



so letting 



2T 



22 18 

Qi Qi 



'3aAi 3a4 
a Ai 



we may write a 2 as 



^B.69) 



3 



°° = ^ no _s - 2a a 2 ) = _ M (B.71) 

Qi 7 Q 2 Qi 7 Q 2 7 



(B.72) 



Ax = 3^ag + ^ag + ^§*3 (B .73) 



(B.74) 



Despite the non-trivial form of the solution we see that it still takes on the 
attractor value at the horizon. 

In terms of the U(l) charges (written implicitly in terms of ao and a 2 ), 
the mass and scalar charge are expressed below 

= 3a 2 - 28a a 3 2 + Z2a\ 

40a aj - 20a 2 a 2 ' 
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2s50aQ a 2 (a — 2a|)(2c4 — 00^2) ^ Qi <3J 

This solution is related to a 3 charge p-brane solution found in [38J - in 
this case we have identified two of the degrees of freedom. 



C Higher Dimensions 



Here we give some more details related to our discussion of the higher di- 
mensional attractor in section 5. The Ricci components calculated from the 
metric, eq. lJHtjj) are, 



r , 2 / n , [d ~ 2)aa 6 „ , ,/-.-,> 

R tt = a la H h aa J (^-l) 

i? rr = -{&(a ,2 + aa") + (d-2)a(a'&' + ao")}/a 2 6 (C.2) 

Rg g = (d - 3) - 2a6a'6' - a 2 ((d - 3)6 /2 + 66") (C.3) 



The Einstein equations from the action eq.f|85jl. take the form, 

(d- 3)(d-3)! a(r) 2 



R tt = uz^nz^ v effWi) ( c - 4 ) 



6( r )2(d-2) a ( r ) 

(C.5) 

where V e // is given by eq. (|HU|) . 

Taking the combination, ^(R rr — ^F-R t t) gives, eq.tjHU). Similarly we have, 

+ a(r) 2 b(r) 2 R rr - (d - 2)R ee 
= -{d - 2){d - 3 - a(r)b'(r)(2a'(r)b(r) + {d - 3)a(r)6'(r))} 
= 2{d r ^fa{r)%{rf - ^ ~ ^ 3) ~ 3) V e// (^) (C.7) 

This gives eq.(jn2J)- Finally the relation, = (d — 3)^?Roe yields, 

(d - 3) 2 (-l + a(r)V(r) 2 ) + 6(r) 2 (a'(r) 2 + a(r)a"(r)) 
+o(r)6(r) ((-8 + 3d)a'(r)b'(r) + (d - 3)a(r)6"(r)) = 0. (C.8) 
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We now discuss solving for the second order perturbation in the metric 
component a, in some more detail. We restrict ourselves to the case of one 
scalar field, 0. The constraint, eq.([92|). to 0(e 2 ) is, 

( d _ 2 )ra' 2 + (d-2)(d- 3)a 2 - 2(0' 1 )V(1 - C-fY^f (C.9) 

-w - w - 3 ) 2 + ^ - v 2 ^M 

H 

+2 (d - 2) ^ m f 2d ThT) {r d H r% + r 3 h r%} = 
r H r 

This is a first order equation for a 2 of the form, 

M + + /a = 0, (CIO) 

where, 

fx = (d-2)r 

h = (d-2)(d-3) 

/ 3 = -2(0' 1 )V(l-(^) d - 3 ) 2 -2(d-2)(d-3) 2 r * " 



r I 1 "'^ r 2(d-3)+l 2 



+2(d-3) 



r 2(d~3) r e-2d 



H 

( A Q ^ ( rf^> /yd y d y 3 \ 

+2(d - 2) [d 6)[r f 2d + r%r%) (C.ll) 



The solution to this equation is given by, 

a 2 (r) = Ce^ - / e~*^dr (C.12) 

where = — J j^dr. It is helpful to note that e T = r{ }_ 3) and, s -^- = . 

Now the first term in eq. (|C.12JL proportional to C, blows up at the hori- 
zon. We will omit some details but it is easy to see that the second term in 
eq. (jC.12|) goes to zero. Thus for a non-singular solution we must set C — 0. 
One can then extract the leading behaviour near the horizon of a 2 from 
eq. (|C.12Jl . however it is slightly more convenient to use eq. ()C.8j) for this pur- 
pose instead. From the behaviour of the scalar perturbation <f>x, and metric 
perturbation, b 2 , in the vicinity of the horizon, as discussed in the section on 
attractors in higher dimensions, it is easy to see that 

a 2 (r) = A 2 (r d ^-r d H 3 ) w (C.13) 
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where, A2 is an appropriately determined constant. Thus we see that the non- 
singular solution in the vicinity of the horizon vanishes like (r — r^)^ 27+2 ^ and 
the double-zero nature of the horizon persists after including back-reaction 
to this order. 

Finally, expanding eq. (jC.12|) near r — > 00 (with C = 0) we get that 
a 2 — > Const + 0(1/ r d ~ 3 ). The value of the constant term is related to the 
coefficient in the linear term for 62 at large r in a manner consistent with 
asymptotic flatness. 

In summary we have established here that the metric perturbation ci2 
vanishes fast enough at the horizon so that the black hole continues to have 
a double- zero horizon, and it goes to a constant at infinity so that the black 
hole continues to be asymptotically flat. 



D More Details on Asymptotic AdS Space 

We begin by considering the asymptotic behaviour at large r of <p\, eq. (|114|) . 
One can show that this is given by 

<M r ) - (^J + (J^J (D.l) 

Here J3/4 stands for a modified Bessel function 15 Asymptotically, I u oc 
r~ 2v . Thus cj) r has two solutions which go asymptotically to a constant and 
as 1/r 3 respectively. 

Next, we consider values of r, r# < r < 00. These are all ordinary points 
of the differential equation eq. ()114|) . Thus the solution we are interested is 
well-behaved at these points. For a differential equation of the form, 

= ^ + P(z)^ + q{z)^ = 0, (D.3) 

all values of z where p(z),q(z) are analytic are ordinary points. About any 
ordinary point the solutions to the equation can be expanded in a power 
series, with a radius of convergence determined by the nearest singular point 

m 

15 Modified Bessel function I v {r),K v (r) does satisfy following differential eq. 

7?l'l{z) + zll(z) (z 2 + v 2 )I(z) = 0. (D.2) 
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We turn now to discussing the solution for a 2 . The constraint eq. (jllO|) 
takes the form, 

2aX + a 2 + (al)'(rb 2 y + ra> 2 = ^/? 2 2 + a 2 ,^^) 2 + ^(r 2 H + ^) + ^ 

(D.4) 

The solution to this equation is given by, 

a 2 (r) = ^ - - / / 3 dr (D.5) 

where 

/ 3 = 2a^' 2 +(ag)'(r6 2 )'+l/3 2 0^ay(9 r ^ 2 -^(^+^)-^ (D.6) 

. We have set the lower limit of integration in the second term at r-g.. We 
want a solution the preserves the double-zero structure of the horizon. This 
means c 2 must be set to zero. 

To find an explicit form for a 2 in the near horizon region it is slightly 
simpler to use the equation, eq. (|109jl . In the near horizon region this can 
easily be solved and we find the solution, 

a 2 oc (r - r H ) 

(27+2) (D7) 



At asymptotic infinity one can use the integral expression, eq. (jD.5|) (with 
c 2 = 0). One finds that f% —> r as r — ► 00. Thus a 2 — > d 2 r. This is consistent 
with the asymptotically AdS geometry. 

In summary we see that that there is an attractor solution to the metric 
equations at second order in which the double-zero nature of the horizon and 
the asymptotically AdS nature of the geometry both persist. 
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